Effective Mass Schrodinger Equation via Point Canonical 

Transformation 

Altu Arda^'G and Rainazan Sever^'IH 

^Department of Physics Education, Hacettepe University, Ankara 06800, Turkey 
^Department of Physics, Middle East Technical University, Ankara 06531, Turkey 

S '. Abstract 

J>-»' Exact solutions of effective radial Schrodinger equation are obtained for some inverse potentials 

by using the point canonical transformation. The energy eigenvalues and the corresponding wave 
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In quantum mechanics, exactly solvable problems are suited in the class having shape 
invariance potential, t^^ The potential in a given class can be converted to one another by using 
the so-called point canonical transformation (PCT) J^^ In the PCT approach, which makes it 
possible to construct a map between the wave equations written for two different potentials, 
it is needed to know the energy eigenvalues, and the corresponding wave functions of a given 
potential (reference potential U{x)) to obtain the energy spectra, and wave functions of the 
other potential (target potential V{x)). The PCT has been used by many authors t^"^' to 
find the energy spectrum, and their wave functions of some potentials in non-relativistic, 
and relativistic scheme. 

Recently, there are some considerable efforts in the literature about the developments 
of the PCT approach to extend to the case of position-dependent mass (PDM) P''^~^^l and 
also the PDM formalism has found a wide application area in the literature. t^^^^^' In this 
Letter, we choose the potentials V{p) = — - + cp^{p = 0, —2) and V{p) = -% + cp^ as the 
reference potential in the case of the radial effective-mass Schrodinger equation (SE), and 
find the exact energy eigenvalues, and their wave functions of the target potentials by using 
the PCT. We follow the approach studied in Ref. [2] where the radial effective mass SE was 
solved for isotropic oscillator. We try to show that the above mentioned formalism could be 
used for inversely linear and inversely quadratic potential forms. 

The SE is written in the case of PDM {niQ = h = I) t^l 

r 1 d"^ 1 dm{x) d \r( \\i ( \ T7 A ( \ n\ 

^2m(x) dx"^ 2m?{x) dx dx ^ " 

and in the case of constant mass 

1 d'^ 
{2 JT - U{y)]i,n{y) = -eni^n{y), (2) 

where En and 0„(a;) are the eigenvalues and the eigenfunctions of the target potential, while 
e„ and ipniy) denote the eigenvalues, and the eigenfunctions of the reference potential. 
By using the transformations f^' 

y = h{x), (3) 

V'n(y) = Tr{x)(j)nix). (4) 

from Eq. (2) we obtain 

{^ + /i(7r(a;);Ma:))^ + /2(7r(x);Ma:)) 

- 2[h'{x)]^ [U{h{x)) - e„] }Mx) = , (5) 



where 



fi{'n-{x);h{x)) = 2tt'{x 

h{n{x)-h{^)) = 71"{X) 



h"{x) 



t:{x) Ti'{x)h'{x) 
1 h"{xy{x) 



(6) 
(7) 



tt[x) -K(x)Tr"(x)h'(x) 

and prime denotes the derivative to the spatial coordinate. If the following equalities are 
satisfied, 



V{x) - E, 



h'(x) = m(x)7r'^{x) , 

[h'{x)r 



m[x] 



U{h{x)) 



f{m{x);h{x)) 



(8) 
(9) 



then the transformation given in Eqs. (3) and (4) is a point canonical transformation, t^l The 
function f{m{x); h{x)) in the above equation is given 



f{m{x);h{x)) 



4:m(x) 



m"{x) 3 I m'{x) 
m{x) 2 y m{x) 



h"'{x) 3fh"{xy 



(10) 



h'{x) 2\h'{x)J 
m{x) and defining the parameter a{x) as a{x) = 4 / h'{x)dx, we obtain 



Taking h'(x) = 

the energy spectrum, the corresponding wave functions, and the target potential as 



Er, 



<Pn{x) 



\ 



h'(x) 



■m(xi 



V'„(/3cr(a;)), 



V{x) = U{/3a{x)) + 



4m(x) 



lm"{x) 7/m'{x) 
2 m(x) 8 y m(x) 



2n 



(11) 
(12) 

(13) 



where the parameter (3 is the scale parameter, t^' 

Giving a PDM function and taking an exactly solvable reference potential U{y), we 
construct a PCT defined by Eqs. (3) and (4) and obtain a target potential, which can be 
used to solve the effective mass SE exactly, eigenvalues, and wave functions of the potential. 

Now, let us look the case of the radial effective-mass SE, which can be written as '^^ 
1 d^ i{i+l) m'{r) (I d' 



L 



\/(r)}0„,(r) = -Ent<pnt{r), (14) 



^2m{r)dr'^ 2m{r)r^ 2m'^(r)\r dr ^ 
where i is the angular-momentum quantum number and the radial wave function is taken 
as ^^rdir) = 0„^(r)/r. 

The radial SE for the case of constant mass reads 



Id^ 
2d^^ 



£'{£' 



U{p)^ilJne'{r) = -ene'ipne'{r). 



(15) 



where £' is the angular-momentum. By using the transformations 



P = h{r), 
ipni'ir) = 7r(r)0„^(r). 



(16) 
(17) 



from Eq. (15) we obtain 



{f-, + Mn{ry,h{r))^ + M7r{ry,h{r))-i'{i' + l)(^ 
^ar'^ ar \ n{r\ 

- 2[h\r)f[u{h{r)) - en,]]i^ne{r) = 0, 



h(r-) 



(18) 



Comparing with Eq. (14) gives 7r(r) = Jh'{r)/m{r) and the following 



V{r) - En, + 



2m(r)r^ m{r) 

m'{r) 



U{h{r)) - tne 



2m{r) \ h{r) 



+ f{m{r); h{r)) . 



2nrfi{r)r ' '' ^ ^ ^' ^ '''' 
which gives the target potential V{r) and its energy eigenvalues Ene for a given radial 
effective mass function m(r)}'^^ 

A) Kratzer potential. 
We consider the following form of the one-dimensional Kratzer potential'^^' 



U{y) 



A B 



(20) 



y y 

which has been extensively used to explain the molecular structure and interactions in 
quantum and molecular chemistry, l^*^' The energy eigenvalues and the corresponding wave 
functions of the reference potential U{y) are'^^' 

2^2 



^ne' = 


2n+l + Vl-165]2' 
= a„yH^^e"«^Fl(-n;l + 2v/=A;2et/);e = 




(21) 


^n£'(r) = 


= V-2en^', 


(22) 



where iFi{—n; g;x) are the hypergeometric type function. Now, we give the explicit form 
of the target potentials and its wave functions for the following three different mass distri- 
butions. 

(1) Mass function m{x) = {6 + x'^y/{l + a;^)^ We obtain a{x) as cr{x) = {1/t)[x + 
{6 — l)/tanx] and the target potential from Eq. (13) 



V{x) 



9(5 



X + 



s-i 

tanx 



1 + 



X + 



S-1 
tanx 



+ 



6-13x^ + 2(2 - S)x^ 
~^ {6 + xY 



(23) 



where 9 = A/ (5 and we have set the potential parameter B as B = A^. The energy spectrum 
and the corresponding wave functions are written from Eqs. (11) and (12) as 

'^ (24) 

^^ exp[-^l3a{x)]iFi{-n;l + 2y/^;2^l3a{x)) . (25) 



En 



16 + x^ 



[X 



x^ 



[(3aix)]- 



(2) Mass function m{x) = a/{6' + x'^) The distribution corresponds to an asymptotically 



vanishing mass function. It gives (7{x) = {\/a/ (5)ln [x + vT+l?]. Therefore, we obtain the 
target potential 

x^ + 25' 



V{x) 



1 



^ ln[x + yl + x^] 
and the energy spectra and their wave functions are 



aln[x + yl + x^] 



la{5' + x' 



2V 



(26) 



(27) 



[X 



[2n+l + Vl-16^'/3']2' 
,^, ^^,y/, [/3a(x)]5+^^exp [-C/3a(x)]iFi(-n; 1 + 2^rA- 2^/3a{x)) . (28) 

where the normalization constant is An = a^'"' a„. 

(3) Mass function m{x) = 1 + tanh{6"x) We obtain the following target potential, 
energy eigenvalues and the corresponding eigenfunctions for this case, 

{6y 7 + tanh{6"x) 



V{x) = e[(T{x)]-^[l + 0[a{x)]-^ 



32 [sinh{6"x) + cosh{6"x)Y 



En 



[X 



[2n + 1 + VI - 16^^/3^J2 ' 
1 + tanh{6"x) [/3a{x)]2 



/^^ 



X exp[-^l3a{x)]iFi{-n; 1 + 2V^; 2^l3a{x)), 



(29) 
(30) 

(31) 

where a{x) = [y/2 / {l35")]tanh^^ [J I + tanh{5"x)/\^] . It is seen that the Schrdinger equa- 
tion can be exactly solved by using point canonical transformation in one-dimension. 

(B) Family of power law plus inverse power potentials. We study the target potential and 
its energy eigenvalues, and wave functions for the following reference potentials in the form 
U{p) = — - + cp^ {p = 0, —2) and U{p) = -^ + cp"^ for the radial SE for the case of position 
dependent mass. 



(1) U{p) 



'- + c. The reference potential has the form 



U{p) 



+ c, 



(32) 



which describes a hydrogen atom with the energy eigenvalue E — c having the energy eigen- 
values and corresponding wave functions, 






c-y(n + £' + l)-^ 



iFi(-n;2£' + 2; 



2a 



n + i' + l 



P)- 



(33) 
(34) 



n + £' + !' 

Taking the radial mass function m{r) = fir'^ and the PCT function h{r) = r" and inserting 
into Eq. (19), we obtain 



V{r) - Ene + 



+ ^)-K~2 <^^'^ y-K-2 ^'^^n£'2u-K-2 



2/x 



/i 



/i 



-K-2 



2/i 



1 



(35) 



- \K'-\{v-l){v-2) + \{v-lf 

which gives two possible solutions 

(a) For the case of z/ = 1 + |; k 7^ —2, from Eq. (35) we obtain 

V{r) = -/iCV-^^t, (36) 

Ene = ~afiC'{n+\C{i)\ + ^-y\ (37) 

where \C{i)\ = =F^^ w4£(£ + 1) + (/t — 1)^ and C is a real coupling parameter. From 
Eq. (17) we obtain corresponding wave functions 



V(,r 



2/i 



-exp 



-(r 



i+# 



.(i+f)(|i:mi+i)+^ 



"V^ + S 
X iFi(-n;2|£(£)| + l;2Cri+t), 



(38) 



where ( = a' /{n + \C{i)\ + |), and 6„ is a normalization constant. 

(b) For the case oi u = k, + 2; k, ^ —2, following the above procedure, we obtain the target 
potential, the eigenvalues and the eigenfunctions by taking into account the independence 
of a potential from any quantum number. 



V{r) 
Ene 



^(2n + 2£(£) + l)', 



(39) 
(40) 



V(r) 



P 



-exp 



(r 



K+2 



,(«+2)(|i:wi+i)-i 



X iFi(-n;2|£(£)| + l;2Cri+2) 



(41) 



where 6„ is a normalization constant and C is the couphng parameter J^' 

The particular case k = —2 is solved by taking h{r) = - ln{ry^^ to obtain the target 
potential. In that case the energy eigenvalues, and the corresponding wave functions only 
exist for the ground state. '^^ We summarize the potential and the energy eigenvalues by using 
Eq. (19) 

(42) 
9i 

- si- («) 

where |A(£)| = ±i^l + 4/iC. The corresponding wave functions are 

0n£(r) = 6„r-i/2[/^(^)]i+|AW| g^^ [-r]ln{r)]iFi{-n; 1 + 2\A{e)\;27]ln{r)), (44) 

where T] = {n+l + |A(^)|)~^ and 6„ = ^(a')"!^^^^! . 

(2) U{p) = —- + -%■ The reference potential has the form 

f/(p) = -i(a + ^), (45) 

and the energy eigenvalues and the corresponding wave functions are listed as '^^1 



V{r) = 

Ent = 


1 rir,-^ r) 


filn{r) 2 



^ni' = - Y (^ + ^ + y (^' + ^)' + 2c) ' , (46) 

V^„^'(P) = anexp[ -^r—rP]iEi{-n;2A + 2;^-^—-p), (47) 

n+A+1 n+A+1 



where A = — | + J(i' + |)^ + 2c. Following the same procedure, we can write the results of 
two possible solutions 

(a) For the case oi u = \ ^ ^]k ^ —2, the target potential and its energy spectrum are 



written as 



\/(r) = -'^r'^-% ^ (48) 

[I 

E - "'^' 



n+^- + ^\CiiW + 2cy\ (49) 



where C{i) = -F^^ y 4^(^ + 1) + (/t — 1)^ — 2c(/t + 2)^ and the corresponding wave functions 
read 



W(r) = hn\—^r''l*exp[ ° r^+2] 

V'^ + S n+i + v/|£(£)|2 + 2c 



X iFi(-n;l + 2v/|£(£)|2 + 2c; ^"^ r^+t). (50) 



where 6„ is a normalization constant. 

(b)For the case oi u = k, + 2; k, ^ —2, we hst the target potential and the energy spectrum 
keeping in mind that any potential function must be independent from any quantum number 



V{r) = -2/iCr'^+^ 



K 



ne 



-2 ' 



where C{i) is given in (a). The corresponding wave functions are written as 

4>neir) = b„ 



(51) 
(52) 



r '^''^exp I 



X iFi(-n;l + 2v/|£(£)P + 2c; 



n+| + y'|£(£)|2 + 2c 
2a 



..K+2\ 



n+i + ./|£(£)|2 + 2c 



(53) 



Finally, we study the particular case where n = —2 by taking h{r) = - ln{r) . By using 
Eq. (19), we summarize the target potential and its energy eigenvalues as 

1 1 ,c C, 



V{r) 



filn{r) l'n?{r) ^ 2 



1 r 1 



^^^ 2 



2 9- 



{--An + - + ^W{^)\+2c\ --} 



(54) 
(55) 



where |A'(£)| = ±|v^l + 4/iC. The corresponding wave functions are 

<Pni{r) = hnV-^'^exp [-r]' ln{r)]iFi{-n; 1 + 2y'|A'(£)|2 + 2c; 2ri'ln{r)), (56) 



where r]' = I n + | + w | A 



'(e\\2 



-1 



2c and 6^ 



a/x. 



(3) f/(p) = -% + cp^. The energy eigenvalues and corresponding wave functions of this 



potential can be written as'^^' 



^nl' 



2n + 1 + -^^(2^ + 1)2 + 8a 



i,ni'{p) = i2cf'p''exp 



:P 



,F,{-n,26 + -;V2^p^). 



(57) 
(58) 



where S = |[1 + \/{2i' + 1)^ + 8a]. In this case, we summarize the two possible solutions as 
(a) For the case oiu = 1 + |; /t 7^ —2, 



V{r) 

Enl 



cv 



/i 



,K+2 



/2cz/4 



/i- 



2n + l + A/2|£(£)|2 + 2a 



(59) 
(60) 



where £(£) = =F^ 
tions 



K+2 



+ 1) + (k — ly — 2a{K + 2)2, and the corresponding wave func- 



(f>neir) = b'^ 



I 2/i 

K + 2 



^5(K+2)+fg^p 



,K+2 



iFi( -n,25 + ^;v^r'^+2)^ ^g^^ 



where 6^ = bn{2cY^'^ and 6„ is a normahzation constant. 

(b) For th( 
eigenvalues as 



(b) For the case oiu = ^ + j, k, ^ —2, we give the target potential and its energy 



V{r) 



-2^Cr" 



2C2/i2c- 



2n+l 



£' 



'(i>\\2 



2a 



(62) 
(63) 



where £'(£) = ^^ 
tions are written as 



+ 1) + A{k, — 1)2 — 2a{K, + 2)2 . The corresponding wave func- 



V(r 



6;jJ/^r^(i+f)+fexj9 



2 



iFi( -n,2(5 + -;v^rtj. (64) 



Now, we proceed the particular case where k = —2 as a last situation by taking h{r) 
- ln{r) . Following the same procedure gives 



V{r) 

Ene 






In^ir) jj, 



I 2c 



2n+l 



+ 2a 



9 
8/I' 



(65) 
(66) 



where T{i) = ±|v^l + SfiQ and O is a real coupling parameter. [21 The wave functions can 
be given as 



?n£ 



(r) = bn^/Jia 



r ^^'"^'exp 



-^^^Mir) 



^F,i -n,26 + l;J^ln'ir)], (67) 



where ,5 = 1[1 + J4 \T{i)\^ + 8a 



In summary, we have used the PCT in the position-dependent mass formalism and choose 
the Kratzer-like potential in one-dimension as the reference potential to construct the exact 
solvable target potential and write the energy eigenvalues, and the corresponding wave 
functions in closed forms. We apply the extended formalism of the PCT to the three- 
dimensional problems and use the potentials V{p) = —- + cp^{p = 0,-2) and V{p) = 
■^ + cp"^ as the reference potentials. The results show that a suitable transformation can 



be used to convert the wave equation into the exact solvable SE and therefore the required 
energy eigenvalues and wave functions are obtained exactly. 
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